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ABSTRACT 

The case is made that the intracluster medium (ICM) in present-day spiral-rich galaxy groups probably has un- 
dergone much slower evolution than that in elliptical-rich groups and clusters. The environments of protoclusters 
and protogroups at z > 2 are likely similar to spiral-rich group environments at lower redshift. Therefore, like the 
ICM in spiral-rich groups today, the ICM in protogroups and protoclusters at z > 2 is predicted to be significantly 
multiphased. The QSO Lya forest in the vicinity of galaxies is an effective probe of the ICM at a wide range of 
redshift. Two recent observations of Lya absorption around galaxies by Adelberger et al. and by Pascarelle et al 
are reconciled, and it is shown that observations support the multiphase ICM scenario. Galaxy redshifts must be 
very accurate for such studies to succeed. This scenario can also explain the lower metallicity and lower hot gas 
fraction in groups. 

Subject headings: galaxies:clusters:general — intergalactic medium — quasars:absorption lines 



1. INTRODUCTION 

Hot X-ray gas (with temperatures of ~ 1 6 — 1 8 keV) is 
found in all galaxy clusters and groups in which it can currently 
be detected (see Burstein & Blumenthal 2002; Mulchaey 2000, 
hereafter MOO). The fact that velocity dispersions of groups 
and clusters correlate well with X-ray temperature has shown 
that galaxy groups are distinct physical entities, not transient 
systems (MOO; Mulchaey et al. 2003). In this paper this hot 
gas is referred to as the intracluster medium (ICM), with the 
understanding that this term also includes galaxy groups in its 
definition. 

Cooling times for hot gas is a trivariate function of its temper- 
ature, metallicity, and density (Sutherland & Dopita 1993). As 
such, in present-day galaxy groups and clusters, hot gas turns 
into cooler, more neutral gas in at least two ways. In the centers 
of clusters, the high density and high metallicity of the hot gas 
combine to predict relatively short cooling times (100's of mil- 
lions of years; Fabian et al. 2002), while in the less dense outer 
parts of clusters, the cooling times are much longer. Separately, 
in groups with low velocity dispersion (i.e., those dominated by 
spiral and irregular galaxies; Mulchaey et al. 1996; Mulchaey 
& Zabludoff 1998; Zabludoff & Mulchaey 1998), the still rel- 
atively metal-rich hot gas can become neutral quickly because 
of its lower temperature and interaction with neutral hydrogen 
in the galaxies in these groups (Burstein & Blumenthal 2002). 
Thus the ICM of spiral-rich galaxy groups today (e.g. the Local 
Group) can be considered to be multiphase. 

The question then is how multiphase is the ICM at high red- 
shift? Answering this question is crucial to understanding the 
evolution of gas in galaxy groups and clusters. In this paper 
the case is made that a multiphase ICM exists in high-z galaxy 
groups and clusters as they are forming (§ 2). The evidence for 
this conjecture that comes from the Lya forest is provided in 
§ 3, where it is shown how to reconcile two recent observations 
that somewhat contradict each other. Further implications of 
the multiphase ICM are discussed in § 4. An = 0.3, A = 0.7 



universe is assumed in this paper, and the distance scales quoted 
below are all comoving. 

2. A LIKELY MULTIPHASE ICM AT HIGH REDSHIFT 

From present-day hot ICM gas, two things must be true when 
groups and clusters are forming: First, because all hot ICM 
gas has a significant abundance of metals (0.1-0.5 times so- 
lar, MOO), much of this hot gas had to have seen the interiors of 
stars. Second, because of the strong correlation of group/cluster 
velocity dispersion with X-ray temperature (MOO), the gas be- 
came hot through heating by supernovae and stirring by the 
virial motion of the galaxies. This observational picture is con- 
sistent with the current hierarchical, clustering, merging sce- 
nario, which predicts that today's galaxies formed from smaller 
objects in an environment dominated by cold dark matter (e.g. 
Navarro, Frenk & White 1995). 

From studies of galaxy groups and clusters in the nearby uni- 
verse (see the Nearby Galaxy Catalog by Tully 1988), the vast 
majority of nearby galaxies are well-associated with either a 
galaxy group or a galaxy cluster. Of the nearly 2450 galaxies 
within 3000 km s" 1 distance of our Local Group that are in the 
Nearby Galaxy Catalog, only 35 are not put in a group, cloud, 
or cluster. And, most of these 35 galaxies are within less than 1 
Mpc of a group or cluster. Furthermore, the majority of galax- 
ies in the Nearby Galaxy Catalog are in galaxy groups (most 
of which are low-velocity dispersion, spiral-rich groups), not in 
elliptical-dominated clusters such as the Virgo cluster. As such, 
it is reasonable to assume that if galaxies are seen at high red- 
shift, they are more likely to be part of a group (even though 
groups are hard to identify at high-z). 

When galaxies are in the beginning throes of formation, the 
ICM should have similar properties to the diffuse ~ 10 4 K in- 
tergalactic medium (IGM), whose temperature scales with the 
density as T oc p a , with a = 0.3-0.6 (Hui & Gnedin 1997). At 
the start of galaxy formation (as well as group and cluster for- 
mation) the density of the ICM may not be much greater than 
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the cosmic mean; hence the ICM cannot be initially much hot- 
ter than 10 4 K. At this stage, the ICM is relatively uniform in 
the gravitational well. 

The ICM will become hotter as galaxies evolve and inject 
hot gas and energy into it. However, the heating before z = 2 
is not strong enough to heat the ICM to 10 6 -10 8 K globally. 
Supernova heating is proportional to the star formation rate, 
which peaks at z — 1-2 (Steidel et al. 1999). Virial shocks 
are not likely to happen at z > 2 or in systems less massive than 
10 n M Q (Birnboim & Dekel 2003). Other sources of heating 
are active galactic nuclei (AGNs), whose number density in- 
creases sharply toward z — 2, but does not evolve much beyond 
Z — 1 (Miyaji, Hasinger & Schmidt 2000). Yet, recent obser- 
vations around protoclusters (Adelberger et al. 2003, hereafter 
A03) imply that AGN heating has not drastically affected the 
ICM at z — 3. Therefore, AGNs may not have a primary role in 
heating the ICM in all clusters and groups. 

On the other hand, cooling of gas at its higher density at 
higher redshift is relatively quick. This is supported by the 
argument that metallicity of this gas had to have been seeded 
quickly by the first generations of star formation. As a result, 
the metallicity of the high-z ICM hot gas probably is not too 
much lower than it is today It is therefore reasonable to assume 
that pockets of moderately-metal-rich, 10 4 — 10 s K gas exist in 
z > 2 protogroups and protoclusters. 

The relatively weak heating and quick cooling necessarily 
lead to a widely distributed, multiphase ICM both in protoclus- 
ters and in protogroups at high redshift. Spiral-rich groups can 
maintain such a distributed, multiphase ICM since they receive 
less energy input from fewer galaxies, and since many of their 
galaxies have neutral gas with which the cooler hot gas can in- 
teract (Burstein & Blumenthal 2002). In other words, the ICM 
in spiral-rich groups today probably has undergone a much 
slower evolution than the ICM in rich clusters and elliptical-rich 
galaxy groups. With this scenario, the environments of proto- 
clusters and protogroups at z > 2 are expected to be similar to 
those in present-day spiral-rich galaxy groups, and therefore the 
ICM in all galaxy systems is multiphase at z > 2. 

In contrast, present-day rich clusters and elliptical-rich 
groups cannot preserve such a distributed multiphase ICM since 
their member galaxies supply large amounts of virial energy to 
the ICM. This picture is consistent with the fact that rich clus- 
ters are mostly observed at z < 1 (Henry et al. 1992), and the 
X-ray emitting ICM in clusters and groups has to be mostly 
formed between z = 2 and today. 

3. THE MULTIPHASE ICM PROBED BY THE LYa FOREST 

The above scenario predicts that foreground galaxies, 
through a strongly multiphased ICM that contains a large num- 
ber of strong Lya absorbers, can influence the Lya forest of 
a background QSO. Accordingly one expects to see a higher 
number of strong Lya absorption lines and a lower mean Lya 
flux at smaller impact parameters between foreground galaxies 
and background QSOs. This gives us a practical criterion: an 
ICM is strongly multiphased if it contains stronger fluctuations 
(measured by Lya absorptions) than the IGM. 

3.1. Observations 

This scenario is consistent with a number of observations. 
Lanzetta et al. (1995), Chen et al. (1998), and Pascarelle et al. 
(2001, hereafter P01) find that galaxies that are not influenced 
by the QSOs (i.e, the "far sample" of galaxies in P01) which are 



at impact parameters p < 180 h~ kpc are much more often as- 
sociated with Lya absorption in the background QSO spectrum 
than those at p > 180 /r'kpc. The rest frame equivalent widths 
W of absorption tend to increase as p decreases from 1 /r'Mpc 
to 10 /r'kpc (P01,Fig. 1). 

It is conservative to assume that most galaxies in the far sam- 
ple are in groups (the same as they are in the local universe). 
The line statistics thus suggest that the ICM in groups does con- 
tain more strong Lya absorbers than the diffuse IGM. The far 
sample in P01 ranges from z = 0.02 to 1, over which groups are 
seen to contain sub-keV gas (MOO; Xue & Wu 2000), consistent 
with a multiphase ICM. On the other hand, Lya lines around 
galaxies that are subject to QSO influence are almost always 
weaker than those in the far sample, and they do not depend 
much on the impact parameter (P01). This shows that the ICM 
cannot remain multiphase under the influence of QSOs, and 
consequently AGNs may be responsible for heating the ICM 
in some groups and clusters at z < 1 . 

At a higher redshift, z — 3, A03 show that mean Lya trans- 
missivity generally increases with increasing impact parameter 
between foreground Lyman-break galaxies (LBGs) and back- 
ground QSOs. It reaches cosmic mean transmissivity of 0.67 at 
p > 6 h~ l Mpc. This result indicates that the ICM at z — 3 has 
a higher overall neutral fraction than the IGM. Since H I con- 
centrates where the Lya lines are, the Lya lines near the LBGs 
must have substentially larger equivalent widths than the Lya 
forest lines. This is consistent with both the line statistics above 
and a multiphase ICM. 

LBGs are more likely to be massive systems (Baugh et al. 
1998) that reside in dense regions, where groups and clusters 
are found today. This strengthens the argument that Lya ab- 
sorption associated with the LBGs out to 6 h Mpc is due 
mostly to the ICM. In fact, the LBGs in A03 are claimed to 
likely be in protoclusters, of which at least a few will evolve 
into rich clusters at z = that contain most of their baryonic 
mass in X-ray emitting hot gas. 

3.2. The Discrepancy 

Despite the general agreement, flux statistics in A03 are not 
consistent with line statistics at p < 0.5 /r'Mpc. An increase 
of mean flux from 0.55 at 1 /r'Mpc to 0.9 at 0.5 /r'Mpc (see 
Figure 1) is reported in A03. 

LBGs are strong UV sources, and the escape fraction of ion- 
izing photons from LBGs at z = 3 is higer than that from lo- 
cal star burst galaxies (Deharveng et al. 2001; Steidel, Pettini 
& Adelberger 2001; Giallongo et al. 2002). It is possible that 
the rise of the mean flux close to LBGs at z = 3 is due to the 
photoionization by the LBGs themselves (similar to the prox- 
imity effect of QSOs, e.g. Bajtlik, Duncan & Ostriker 1988; 
Giallongo et al. 1996; Scott et al. 2000), which is absent in 
low-redshift galaxies. One can roughly estimate the size of 
influence by LBGs as follows. If the luminosity of a typical 
QSO is 10 46 erg s"' and if it has a line-of-sight (LOS) prox- 
imity zone of 36 h~ l Mpc (4000 km s" 1 at z = 3), then an LBG 
with a luminosity of 10 41 erg s~' would create a highly ion- 
ized (higher than the IGM) bubble with a radius of at most 0.8 
/r'Mpc, where a proportionality between volume and luminos- 
ity has been assumed. The radiation of QSOs is most likely 
anisotropic, as suggested by null detections of foreground QSO 
proximity effect (Crotts 1989; M0ller & Ksergaard 1992; Liske 
& Williger 2001), so that the relatively isotropic radiation from 
LBGs should considerablly reduce its radius of influence. In 
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addition, QSOs have harder spectra than LBGs, which help to 
create larger proximity zones, owing to larger fractions of high- 
engergy photons that can travel farther to ionize H I because 
of their smaller cross sections for photoionization. Therefore, 
an estimate of 0.8 /r'Mpc is a conservative upper bound for 
the radius that can be affected by radiation from LBGs. Al- 
ternatively, A03 suggests that supernova-driven winds with a 
sustaining speed of 600 km s" 1 for a few hundred million years 
may be able to deplete H I and reduce Lya absorption in a re- 
gion of comparable size. 

Regardless how the highly ionized bubble is created, how- 
ever, Weinberg et al. (2003) show with hydrodynamical simu- 
lations that even if one removes all H I within 1.5 /r'Mpc from 
LBGs in real space, the mean flux at p < 0.5 /r'Mpc is still 
lower than 0.67. It is because the infalling H I from large radii 
could be easily seen within 0.5 /r'Mpc in redshift space. This 
is similar to finger of god effect, which smears real-space in- 
formation along the LOS. For example, the Hubble constant at 
z = 3 is 446/z km s" 1 Mpc - ' , so that an infall speed of 1 20 km s" 1 
will enable Lya absorbers at a LOS distance of 1.5 /r'Mpc 
from an LBG in real space to appear less than 0.5 /r'Mpc from 
the LBG in redshift space. Since the peculiar velocity of a gas 
clump in a cluster easily exceeds 120 km s~', it will completely 
erase the LBG bubble in redshift space. 

In summary, the flux statistics found by by A03 are incom- 
patable with the line statistics found by P01 at p < 0.5 /r'Mpc. 
This difference cannot be easily attributed to ionization by 
LBG themselves, through either radiation or winds, unless 
LBGs could ionize an unrealistically large volume (radius » 
1 /r'Mpc) that cannot be erased by the virial motions of gas 
clouds. 

3.3. Resolving the Discrepancy: Redshift Errors 

The mean flux at an impact parameter from a galaxy is mean- 
ingful only if it is calculated within a compatible distance from 
the galaxy along the LOS. The rms error of galaxy redshift is 
a z = 0.002 in A03, which corresponds to 1.3 /r'Mpc at z — 3. 
In turn, redshift errors can lead to an inaccurate estimate of the 
mean flux within p = 1.3 /r' Mpc. A03 evaluates the mean 
flux in their data using the flux-galaxy correlation function to 
suppress the effect of redshift error. This method only helps 
if the sample is large. However, A03 has only 3 galaxies at 
p < 0.5 /r'Mpc without Ly/3 contamination (3 more with Ly/3 
contamination), so redshift errors could still alter mean flux 
within p = 1.3 /r'Mpc. 

The effect of a small redshift error Az on the mean flux 
is examined using the Keck high resolution spectrum of QSO 
HS 1700+64 (for observational details and data reduction tech- 
niques, see Kirkman & Tytler 1997). This QSO is at z = 2.74, 
and the spectral resolution of the data used is 8 km s"' . Only 
the segment between z = 2.16 and 2.70 is analyzed in order to 
exclude Ly/3 absorptions and the QSO proximity zone. In an 
idealized case that a galaxy is at the center of a Lya line, mean 
flux of the Lya forest around the galaxy redshift should increase 
with the scale over which this mean flux is calculated, until it 
reaches the cosmic mean. 

A mock sample of 59 galaxies is made by assigning the po- 
sitions of strong Lya lines (4 A > W > 0.5 A) to galaxies. If 
two lines are closer than 2.5 A, the stronger one remains. The 
line width criterion is arbitrary, but it only contributes to the 
dispersions and details, not the general trend. Line widths are 
adopted from Dobrzycki & Bechtold (1996) with corrections to 



line positions. Impact parameters are not assigned; rather mean 
fluxes are calculated within the LOS distances (still reported as 
p) that are equal to the impact parameters as argued above. 

Mean flux around accurate galaxy positions is shown in Fig- 
ure 1 with open circles, which decreases monotonically as p 
decreases. A small systematic error of Az = 0.001 (filled cir- 
cles) can easily raise the mean flux at p < 1.5 /r'Mpc. If the 
error is Gaussian with a dispersion a z = 0.002 (solid line), the 
mean flux drops slightly from p = 2 /r'Mpc to p = 0.5 /r'Mpc, 
but it never rises at smaller p. 

The same Gaussian error was tested by A03, but no signif- 
icant effect was found. If the positions of the few galaxies in 
A03 at p < 0.5 /r'Mpc are all off by Az = 0.002, one cannot 
recover a meaningful ensemble average of the mean flux by 
adding the Gaussian error <r z = 0.002 to their positions. The po- 
sitional error of a small number of galaxies is simply not well- 
characterized by a Gaussian distribution. More observations are 
needed to conclusively establish the mean flux at small impact 
parameters. 

If the deficit of Lya absorption at p < 0.5 /r'Mpc is the re- 
sult of errors in galaxy positions, then the flux statistics and the 
line counting statistics give a consistent picture of a multiphase 
ICM that exists in high-z galaxy agregations and low-z groups. 
If not, one will have to explain how LBGs ionize a sphere of 
radius 3> 1 /r'Mpc in real space. 

3.4. Neutral Hydrogen in the ICM at z — 3 

Given the mean flux, (F), one can estimate the density of 
neutral hydrogen n H using (Peebles 1993) 

n H ~2.4x 10""f7' /2 /z(l+z) 3/2 ln( J F}-'. (1) 

The H I column density Nh at z = 3 is roughly 

N H ?t n H l ~ 8.1 x lO'Mn^}"' (///r'Mpc) cm -2 , (2) 

where / is the (comoving) LOS path length through the ICM. 
With the lowest mean flux, (F) = 0.52, near LBGs in A03, and 
/ = 12 /r'Mpc for a spherical ICM that extends to 6 /r'Mpc 
(§3.1), equation (1) yields an upper limit of Nh = 6.4 x 10 14 
cm" 2 in the ICM at z = 3. This is several orders of magnitude 
lower than Galactic values, so it can easily escape detection 
from 21 cm radio and other observations, even if it still exists 
in the ICM today. 

4. DISCUSSION AND CONCLUSIONS 

The Lya forest in the vicinity of foreground galaxies is po- 
tentially a great probe for studying the ICM. Caution must be 
taken, however, since knowing the accuracy of galaxy redshift 
is crucial when probing at small impact parameters. Never- 
theless, the available observations point to a multiphase ICM at 
z = 3, which exists partially in the form of strong Lya absorbers. 
It is also evident from the Lya forest that a multiphase ICM 
is preserved in spiral-rich galaxy groups today, consistent with 
what can be predicted about the cooling times for the hot gas 
in such galaxy groups (Burstein & Blumenthal 2002). An inde- 
pendent piece of evidence comes from recent discovery of the 
warm-hot 10 5 -10 7 K intiagroup medium through observations 
of the high velocity O VI absorption clouds around our Galaxy 
at a likely distance > 100 kpc (e.g., Nicastro et al. 2003). For 
rich clusters, the ICM may have been homogenized by galax- 
ies, so that it is unlikely to be multiphase except in core regions 
where cooling times are short, and therefore expected to show 
much less Lya absorption than the IGM. 
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It is found that ionization by LBG themselves via radiation 
or supernovae-driven winds cannot easily account for the lack 
of Lya absorption at p < 0.5 /r'Mpc in A03, because LBGs do 
not have sufficient energy to ionize a region of size » 1 /r'Mpc 
in real space. Errors in galaxy positions can cause the apparent 
rise of mean flux at p < 0.5 /z~'Mpc found by A03. 

The origin of the multiphase ICM is the competition between 
heating and cooling in these environments. The detailed prop- 
erties of the absorbers are not well determined by Lya obser- 
vations alone. This is because even though cool (~ 10 4 K) and 
dense absorbers are common at z = 3, shocked gas clumps ap- 
pear more often at low redshift. These gas clumps can have 
high H I column densities (e.g. Dave et al. 1999), in spite of 
higher temperatures (> 10 5 K). 

The existence of a multiphase ICM requires effective cooling 
and low heating. This is consistent with groups having lower 
fraction of mass in hot gas than clusters, if a higher fraction of 



the ICM in groups is kept in cool and warm phases. In addi- 
tion, since galaxies heat the ICM by feedback and motion, less 
heating means lower metallicity in the ICM in galaxy groups, 
which is observed (MOO). 

A multiphase ICM may also contribute to the excess of en- 
tropy in the cores of poor clusters and groups (Ponman, Cannon 
& Navarro 1999), which requires lower electron density than 
predicted by the hierarchical clustering model. Efficient cool- 
ing in poor clusters and groups may reduce hot gas density in 
a short time. If this picture is correct, one might circumvent 
the problem of missing soft X-ray emission encountered by the 
cooling flow model (Fabian et al. 2001; Wu, Fabian & Nulsen 
2001), which lowers the hot gas density by on-going cooling. 

We thank D. Tytler for kindly providing the Keck spectrum 
HS 1700+64, and the referee for helpful comments. 
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FIG. 1. — Mean Lya flux around galaxies. The open circles, filled circles, and the solid line are the cases with no error in redshift, a systematic error of Az = 0.001, 
and a Gaussian error with dispersion a, = 0.002, respectively. The error bars of the three are similar, but only plotted for the open circles for clarity. The diamonds 
are from A03. The segment of the Lya forest in our analysis is from z = 2. 16 to 2.70, so the cosmic mean flux is higher than that at z = 3. This leads to the difference 
between A03 and the rest of the data at large impact parameters. 



